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Competing scenarios for quantum critical points (QCPs) of strongly interacting Fermi systems 
signaled by a divergent density of states at zero temperature are contrasted. The conventional 
scenario, which enlists critical fluctuations of a collective mode and attributes the divergence to 
a coincident vanishing of the quasiparticle strength z, is shown to be incompatible with identities 
arising from conservation laws prevailing in the fermionic medium. An alternative scenario, in which 
the topology of the Fermi surface is altered at the QCP, is found to explain the non- Fermi- liquid 
thermodynamic behavior observed experimentally in Yb-based compounds close to the QCP. It is 
suggested that combination of the topological scenario with the theory of quantum phase transitions 
will provide a proper foundation for analysis of the extended QCP region. 
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a. Introduction. Fundamental understanding of the 
behavior of Fermi systems in the vicinity of quantum 
phase transitions persists as one of the most challenging 
objectives of condensed-matter research. As it involves 
second-order transitions occurring at a critical density 
Pc^ the problem is even more difhcult than in the clas- 
sical regime, since the description of quantum fluctua- 
tions entails a new critical index, the dynamical critical 
exponentjii^ In several heavy-fermion metals — notably 
Yb-based compounds^ — critical temperatures Tn{B) can 
be driven to zero by extremely weak magnetic fields B, 
creating a quantum critical point (QCP). It is commonly 
believed that low-temperature fluctuation contributions 
to the free energy, specific heat C{T), and other ther- 
modynamic quantities must then follow power laws in T, 
the Somnierfeld ratio C{T)/T being divergent at T — > 0. 

From a pedestrian standpoint, such non- Fermi-liquid 
(NFL) behavior must extend some distance from the 
QCP, implying separation at T = of a domain of mag- 
netic ordering from a Fermi-liquid (FL) regime — as is in- 
deed the case in the heavy-fermion metal YbAgGci^ How- 
ever, this example is unique; otherwise, the two phases 
seem to abut each other at the quantum critical point 
(QCP). Since the standard FL formalism is applicable 
on the FL side of the QCP, C{T)/T is proportional to 
the effective mass M* in this region and it follows that 
M* must diverge at the QCP. 

Conventional arguments that quasiparticles in Fermi 
liquids "get heavy and die"^ at the QCP commonly em- 
ploy the textbook formula 
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where 



ppjM and the derivative is evaluated at 



p = Pf and e = 0, single-particle (sp) energies be- 
ing referred to the chemical potential /z. The factor 
z = [1 — (9S(p, £)/(9£)q]~^ is the quasiparticle weight of 



the sp state at the Fermi surface. The conventional be- 
lief, traced back to Ref. 0, holds that the divergence of 
M* at the QCP is caused by the vanishing of the z factor 
in Eq. ([T]), stemming from the divergence of the deriva- 
tive {dYi(jp,£] pc) I de)f^ at implicated second-order phase 
transition points. 

However, this scenario is problematic. As will be seen, 
the z-factor does not vanish at the points of second-order 
phase transitions. It will be argued that the divergence 
of the density of states N{T) at the QCP is in fact asso- 
ciated with a rearrangement of single-particle degrees of 
freedom^ rather than with critical fluctuations. Even so, 
the divergence of N{T) at the QCP does give rise to some 
second-order phase transitions, occurring at T = Tn in 
the vicinity of the QCP. Accordingly, the full pattern of 
phenomena in the temperature interval from to T > Tn 
in the QCP region is determined by an intricate inter- 
play between the two mechanisms. The present analysis 
is limited to the disordered side of the QCP where the 
impact of sp degrees of freedom is decisive, while the role 
of critical fluctuations is suppressed. Other regions of the 
phase diagram will be analyzed elsewhere. 

b. Fault lines of the conventional scenario for the 
divergence of M* . We begin by exposing inconsisten- 
cies of the standard derivation leading to divergence of 
(9S(p, e) /9e)g in the vicinity of second-order phase tran- 
sitions. This derivation is based on a fundamental rela- 
tion of many-body theory. 
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where U is the totality of diagrams of the scattering am- 
plitude, irreducible in the longitudinal particle-hole chan- 
nel. 

The treatment in question retains only the pole part 
G"' of the sp Green function G — zC^ + G^' and a singular 
part of the diagram block U that is supposedly responsi- 
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ble for the divergence of {dT,{p, e)/de)Q. By definition of 
the block U, there is no the direct spin-fluctuation contri- 
bution to e) /9e)g , while the exchange term has the 
form UaS-ff3{p, sAi^i) = g'^(^a/3cr-fSX{<l,<lo) involving the 
spin susceptibility go)j where q — \p — l\, qo = e — ei, 
and 5 is a dimensionless effective coupling constant. 

On the FL side of the QCP, Tias{p,£) — ^{p,£)6aS 
and ImG"^(p, e) = —2Trsgn[e{p)]S{e — e{p))- In accord 
with the FL perspective and conventions, we now write 
UaS-yp = Uo5a5S.^f3+Us(TaS<T-ip, and the single component 
Uo enters Eq. yielding 



de 



9G(/,£i) d\dei 
dei (27r)4 



(3) 



Applying the identity (Tapcr^s = ^SasS^p - ^CTaStr^p, we 
have Uo{q, qo) = {3/2)g'^xiq, Qo) and Eq. Q reduces tc^-^ 



\ de 



dp \ f f i.-,<ldq 



(4) 

A key assumption made in Ref. l6j, and generally 
adopted in subsequent treatments, is the Ornstein- 
Zernike (OZ) form 

x{q, go = 0) = xoz(g) = , n (5) 

C -I- g 

for the static correlation function xil)^ with the corre- 
lation length ^ diverging at the critical point. Inserting 
this ansatz into Eq. (g]) together with vq = {de{p)/ dp)^ = 
VpM/M*, one arrives at 
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(6) 



which implies that {dll{p,e] p) / de)^ diverges at ^ ^ oo. 
However, the applicability of the OZ approximation to 
evaluation of £)/9£)g in the QCP domain has 

never been proved. 

This deficiency exhorts us to check the compatibility 
of the OZ approximation in homogeneous matter with a 
set of identities involving the derivative 
all having the same structure as Eq. ([2]) • In so doing, we 
observe that Eq. ([2]), which follows from particle- number 
conservation with the aid of the scalar gauge transfor- 
mation '^{t) — > ^'(t)e*^*, is but one instance of a class 
of similar identitiesi^ Any conservation law existing in 
the medium generates a corresponding identity involving 
(9S](p, e)/9e)g. For example, momentum conservation in 
homogeneous, isotropic matter, associated with the vec- 
tor gauge transformation — > \E'(i)e'P'^', results in 
the well-known Pitaevskii relatior 
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Uo[p,eXei) 



dG[l,e i) (pi) dldei 
dei 



(7) 



p2 {2t:)H 

Analogously, in the model of Ref. 6 the spin operator <t^ 
commutes with the Hamiltonian, and the gauge transfor- 
mation vl/(t)e*'^3^* leads to the relation 

(9S(p,e) r , , ,5G(/,ei) d\dei 



de 



Us{v,e,\,ei)- 



dei (27r) 



(8) 



Even more conservation laws exist in nuclear and dense 
quark matter, each providing an identity like Eq. ^ . 

The standard manipulations applied to relation ([2]), 
leading to the result ^ via ansatz ([5]), can now be re- 
peated for any such conservation identity. Irrespective of 
which identity is chosen, a divergent result is obtained 
for {d'E(p, e; pc)/de)Q. Importantly, however, the signs of 
the divergent components of this derivative do depend on 
the choice made. For example, in the case of critical spin 
fluctuations, adoption of the anzatz ([5]) results in diver- 
gence of the derivative {dT,{p, e; Pc)/de)g whether Eq. ^ 
or Eq. ([5]) is adopted, but different signs are delivered, 
since in the OZ approximation the blocks Uo and have 
opposite signs. Thus, Eq. ([3|) provides an "acceptable" 
negative sign, whereas Eq. fH]) gives a "wild" positive sign 
(and a meaningless limit for z). In the case of short- wave- 
length fluctuations with critical wave number qc, a similar 
wild result is obtained from Eq. ^ because the prefactor 
of the divergent part of {d'S{p, e)/de)Q differs from that 
derivedii from Eq. ([3]) by a factor cos 0c = 1 — Icf^P^- 
Since the nonsingular components of the block U are in- 
capable of compensating the divergent OZ contributions 
to {dY^{p,e; pc)/de)Q, we must conclude that the result 
([6]) is fallacious, and that more sophisticated methods 
must be applied to clarify the situation in the critical- 
point region. 

We call attention here to the situation for classical 
second-order phase transitions, where the OZ correla- 
tion function (O is altered by scattering of the fluc- 
tuations themselves. 12, As a result, the actual correla- 
tion function ^(^j Pc) decays more rapidly at large dis- 
tance r than xoz [f ^ Pc) cx: l/r. In momentum space, 
x{q 0,gc) behaves^- as 1/q'^^'^, with 77 > 0, com- 
pared with xoz{q) cc 1/q^. If a similar alteration of 
x{q — * 0, Pc) occurs at r ^ 0, then the integration lead- 
ing to ([6]) is saturated at p pc, ensuring that z{pc) ^ 0. 

c. Topological scenario for the QCP. With the con- 
dition z{pc) = ruled out, the effective mass in Eq. ([T]) 
can only diverge at a density poo where the factor 
in square brackets, or equivalently the group velocity, 
changes sign. Such a QCP can be examined based on 
the FL equatiortiS 



de{p) del f ,5n(pi) , 



(9) 



where ep is the bare sp spectrum and dr is the volume 
element in 3D or 2D momentum space. The T = group 
velocity, being a continuous function of the interaction 
function /(p,pi), changes its sign on the Fermi surface 
at the critical density poo- In 3D homogeneous matter 
the critical condition is 



1 = /l(Poo,Poo)PooM/37r^, 



(10) 



where /i is the first harmonic of / and Poo — (37r^poo)^/^. 
In this scenario, the QCP is associated with a rearrange- 
ment of single-particle degrees of freedom; no collective 
parameter is involved, and the symmetry of the ground 
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state is not broken. Such topological phase transitions, 
induced by the interactions between quasiparticles, have 
been discussed for over two decades.— 

For a homogeneous medium there are in general two 
ways to reahze a divergent density of states iV(0,poo)- 
Both options are associated with bifurcation points pb 
of the equation e{p, poo) = 0. As a condition for the 
divergence of the effective mass M'*, Eq. (jlOp refers to 
the caseiS p,^ = pp in which the sp spectrum e{p) has an 
inflection point. In the second option, where pb ^ Pfi M* 
remains finite, while N(0,pcc) diverges due to vanishing 
of the group velocity at the bifurcation point. 

Thus far we have dealt only with homogeneous sys- 
tems. A comparable analysis of NFL behavior of heavy- 
fermion metals must include the effects of anisotropy, 
which are of special importance in the QCP region. An 
early study of topological phase transitions in anisotropic 
electron systems in solids, induced by electron-electron 
interactions, was carried out in Ref. il4l . 

Here it will be instructive to address the 2D electron 
liquid in a quadratic lattice, assuming the QCP electron 
Fermi line to be approximately a circle of radius Poo, 
with the origin shifted to (7r/a,7r/a). Since the group 
velocity w„(p) = de{p)/dpn now has a well pronounced 
angular dependence, the topological anisotropic QCP is 
to be specified by the vanishing of f„(p, 0, T = 0, poo) at 
the single point p = Poo,4> = 0- On the FL side of the 
QCP where dvn{p, <j>)/d<l) > 0, one has 

Vnip, (f>;T = 0,p) = Vn{p,(j) = 0) + a^cj)^ + ap{p - poc), 

(11) 

with Vnip, (f> = 0) = anip — Poo)^ as in the inflection-point 
case treated in Ref. Hj. The QCP density of states 



iV(T,p)cx / n{e)il-n{e)) 



de ( 



vn(p{^)A;T, p)' 



(12) 



which determines the specific heat C{T) — TdS/dT ^ 
TN{T) and thermal expansion coefficient P{T) ^ 

~dS{T,p)/dP TdN{T,p)/dp, is evaluated utilizing 

Eq. (|11|) . Following analytic integration over 0, we have 



C{T,p^) (X / n(e)(l-n(e))- 



de 



v'J\p{e),<t> = Q-T) 



(13) 



/3(T,poo)cx / n(e)(l-n(e))^— ^^^^^22^^ « 0(1), 

(14) 

where v'„{p) = dvn{p)/dp — ap by Eq. pT|) . In arriving 
at the overall temperature behavior, we have introduced 
the estimatei^ Vn{p,(t> — 0',T) oc T^/"^ stemming from 
Eq. pT|) and valid in the relevant interval e ~ T. The 
results ([T^ and may be combined to determine the 
behavior of the Griineisen ratio, 

r(T, Poo) = /3(T, Poo)/C(T, Poo) cx (15) 

at variance with the FL result F oc 0(1). 



Imposition of an external magnetic field greatly en- 
larges the scope of challenging NFL behavior, as re- 
flected in the magnetic Griineisen ratio T^iig{T, B) = 
~{dS{T,B)/dB)/C{T,B). Wc now analyze this key 
quantity within the topological scenario, again follow- 
ing the path established in Refs. ITslflsl . The original 
Fermi line is split into two, with consequent modifica- 
tion of field-free relations such as (O and (fT^ through 
the appearance of the sum of quasiparticle occupancies 
n±{e) = {1 -I- exp [(e ± pei3)/T]}~"^. As a consequence, 
v^^ is replaced by the sum of quantities {de{p, 4>) / dpn) ^ 
evaluated at e(p, (f>) ± fieB. Analytic integration over (jj 
still goes through and yields the sum of square roots of 
these quantities. In the limit T ^ PeB, terms in this 
sum linear in r = p^B/T cancel each other, such that 
the net result is proportional to , leading to 



rmag(7' > PeB) OC T 



(16) 



More specifically, upon integration over (j) in the field- 
perturbed formulas for C{T,B) and S{T,B), the ensu- 
ing expressions involve the sum l/{e + p^B)^^/^ + l/(e — 
PeB)~^^^ , multiplied by a factor containing only n(e). In- 
tegrating over the dimensionless variable y = e/T, the ra- 
tio S{B,T)/C{T) is determined as a function of only, 
yielding the result (|16p at r ^ 1. 

In the opposite limit r ^ 1, the density of states 
N{T = 0,B) diverges at a critical magnetic field -Boo, 
where the function Vn{p,4i,T = 0, -Boo) vanishes on one 
of the two Fermi lines p^ {(f) specified by 



e(p±,0)±MeB = O. 



(17) 



The field-induced splitting that alters the relevant Fermi- 
surface group velocity can be compensated — for example, 
by doping — thereby providing the means for driving the 
system toward the QCP. 

At -B > -Boo, the critical quantity Vn{p,(t> — 0;T — 
0, -B) = Vn{B) becomes positive and FL behavior is re- 
covered, as in the isotropic caseJ^ii^ To evaluate the 
critical index specifying the divergence of the density of 
states N{T = 0, -B ^ -Boo) oc Vn^^'^(B), we calculate the 
spectrum from Eq. ([T7]) [as in Ref. [lj| and insert the re- 
sult into Eq. (fT7|) . obtaining Vn {B) oc{B- B^fl^ and 
N{T = 0, -B) oc (B - B^)-^/^. Thus 



C{T -^0,B) = S{T ^ 0, B) oc T{B - B^ 



-1/3 



(18) 



and (3{T ^ 0) oc T{B - Boo)^^, so that r(T ^ 0) oc 
(-B — Boo)^^^^- Importantly, we arrive at 



Tm.g{T ^ 0) ^ ^{B - B^y 



(19) 



Such a divergence was first predicted within scaling 
theory,. ^-fi, in which the peak of r„jag{T = 0, B) is located 
at Be, the end point of the line T]y{]3) where T/v(-Bc) = 0. 
In the topological scenario, -Boo may or may not coincide 
with Be', this issue will be pursued in further work. 
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d. Discussion. The results (IT3l)-((T6l), and 
(|19p are in agreement with available experimental 
dataji^iiiii^'i^ obtained by the Steglich group in com- 
prehensive studies of the thermodynamic properties of 
Yb-based heavy-fermion metals. These data also pro- 
vide a test of modern phenomcnological scaling theo- 
ries of the QCP i^°i^^ The outcome of this test, as aired 
in Refs. ,18:,2l], is that no single model based on 2D or 
3D fluctuations can describe these data, which require 
the following set of critical indexes having low probabil- 
ity: dimensionality d — I, correlation- length exponent 
v = 2/3, and dynamical exponent 7 = 3/2. 

Recent studies of peaks in the specific heat C{T, B) 
in Yb-based compounds reveal another difficulty con- 
fronting the phenomcnological theory of second-order 
phase transitions in the QCP region. According to this 
theory, at B = the difference T — Tm is the single 
relevant parameter determining the structure of the fluc- 
tuation peak of the Sommerfeld ratio C{T)/T. How- 
ever, comparative analysis of corresponding experimental 
d&iB^^ in YbRhaSia and YbRh2(Sio.95Geo.o5)2 shows 
clearly that the structure of this peak is not universal. 
As the QCP is approached, the peak gradually shrinks 
to naught, again bringing into question the applicability 
of the spin-fluctuation scenario in its vicinity. 

Thus, while the spin- fluctuation mechanism remains 
applicable at finite T ~ Tn{B), it becomes inadequate at 
the QCP itself. Accordingly, the relevant critical indexes 
of scaling theory must be inferred anew from appropri- 
ate experimental data on the fluctuation peak located at 
Tn(B). Furthermore, the existing description of ther- 
modynamic phenomena in the extended QCP region, in- 
cluding the peak at Tn{B), must be revised by integrat- 
ing the topological scenario with the theory of quantum 
phase transitions,^ 

The posited suppression of critical fluctuations in 
YbRh2(Sio.95Geo.o5)2, which is situated extremely close 
to the QCP in the sense that Be = 0.027 T, conflicts with 
the conclusion of Ref. [13 that the system is on the verge 



of a ferromagnetic instability. This assertion is based on 
extraction of the Stoner factor from measurements of the 
Sommerfeld- Wilson (SW) ratio i?sw(r) oc x{T)/C{T). 
However, such an extraction is straightforward only in 
homogeneous matter, where the magnetic part of the 
Hamiltonian is specified by the Bohr magneton /i^. In 
dealing with electron systems of solids, this strategy is 
inconclusive unless a reliable replacement fj,es for /i^ is 
known. The authors of Ref. [23 have chosen the effec- 
tive Bohr magneton fi^e to be 1.4-1.6 fiB, as determined 
from data on the magnetic susceptibility itself. Such a 
choice suffers from double counting. If instead one uses 
the value ^^s = 4.54/iB appropriate for the atomic state 
of Yb^+, then the Stoner factor derived from the data 
remains below 3, and the conflict is resolved. 

e. Conclusion. The conventional view of quantum 
critical phenomena, in which the quasiparticle weight z 
vanishes at points of related T — second-order phase 
transitions, is incompatible with a set of identities based 
on gauge transformations associated with prevailing con- 
servation laws. We have traced the failure of the standard 
scenario to the inapplicability of the Ornstein-Zernike 
form x~^{(l) = + C '^ for the static correlation func- 
tion x('J') in the limit ^ — > 00. We have discussed an 
alternative topological scenario and demonstrated that 
its predictions for the thermodynamics of systems on the 
disordered side of the QCP are in agreement with avail- 
able experimental data. Based on these data, we infer 
that close to the QCP the role of single-particle degrees 
of freedom is paramount, while the effects of critical fluc- 
tuations build up on the ordered side as the system moves 
away from the QCP. 
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